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Abstract 

o 

£SJ . In this paper we study Hawking-like radiation in a FRW- Universe using semi-classical tunnelling 

approach and the Hamilton- Jacobi method. Radial null geodesies are used to picture the process 
as a "tunnelling" of particles from behind the apparent horizon and Hawking -like temperature 
has been calculated. Quantum corrections have been evaluated in the Hamilton- Jacobi Method by 
solving Klein-Gordan wave equation and the temperature agrees at the semiclassical level. Also it 
CNI ' is found that Hamilton- Jacobi formalism does not depend on the choice of the coordinate system. 

Finally, leading order corrections to entropy has been calculated. 
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. 1 Introduction 

CO 

In recent years there are lots of works [TJH1I31I11I5JIE1II1IEJIH1 HQ1 E] m universal thermodynamics with a view to 
\& • establish the equivalence between the Friedmann equations and the first law of thermodynamics. These studies 

are mostly concentrated to FRW model of the Universe (For universal thermodynamics in LTB model see ref 
12 ). In these works universe bounded by the apparent horizon is chosen as a thermodynamical object and the 
entropy and temperature of the boundary are assumed to be 



5=A T =^— 
4G ' 2irR A 

where A and Ra are the area and radius of the apparent horizon respectively. Jacobson [13] first formulated 
the Einstein field equations using the fundamental Clausius relation dQ = TdS. Subsequently, recent progress 
of black hole thermodynamics has been applied to the universal thermodynamics [TU [T5J IE1 El [TB] . 

To study the Hawking radiation there are two basic approaches — Firstly, the tunnelling approach of Parikh 
and Wilczek 021(201 HI] and secondly the standard Hamilton-Jacobi (HJ) method (also known as complex path 
integral formalism) [22, 23:. Due to simplicity of its calculation there are lot of works [24] [25J [26J [27l [28] in this 
method. Recently these approaches has been applied to universal thermodynamics. Hawking-like temperature 
is calculated using both the methods at the apparant horizon of the FRW universe [3H1 1301 1311 I3H 1331 131] (see 
also ref [35]) for LTB model of the universe). Initially there was a problem of factor 'two' [311 1321 1331 134] in the 
Hawking-like temperature and it was resolved subsequently considering contribution to the probability from the 
time co-ordinate upon crossing the horizon j36j [371 131] ■ 

In the present work, we make an exhaustive study of the FRW universal thermodynamics using both the 
methods. We review the radial null geodesic method using Kodama vector in section 2. Section 3 deals with 
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HJ approach in which full quantum corrections have been evaluated. Semiclassical Hawking-like temperature 
is also calculated using Painleve type coordinates by HJ formulation in section 4 and the temperature agrees 
with the previous results. Finally, the paper ends with a brief conclusion in section 5. 



2 Hawking-like temperature by radial null geodesic method:An Overview 

In recent years the popular approach of evaluating Hawking radiation is the radial null geodesic method 
[TpJ HDJ HI] (known as semi classical tunnelling analysis). The method is simple compared to HJ method but it 
has some limitations, namely 

(i) The method is applicable only for massless particles 

(ii) one has to use only Painleve type coordinates to avoid singularity at the horizon 
(Hi) there is a discrepancy of factor two in this method 

(iv) there is no general method to include quantum effects (back reaction). 

The basic idea in this semiclassical approximation (WKB approximation) is that the emission rate for the 
s-wave emission of a massless particle can be related to the imaginary part of the action of a system. In the 
s-wave approximation, particles can be viewed as a massless shell moving along radial null geodesic. Now 
compared to static black hole cases, there is one basic difference as in the present case the metric coefficients 
depend on both radius and time . Hence there is no longer time translation Killing vector field. We shall have 
to use the Kodama vector |39j which is time like inside the horizon and the associate energy. 

The homogeneous and isotropic model of the universe described by the FRW metric as 

ds 2 = -dt 2 + ^-^dr 2 + a 2 (t)r 2 dn 2 2 (1) 
1 - kr 2 w * 

Here the FRW coordinate (t, r, 9, cf>) is an orthogonal co-moving coordinate with 't' the co-moving time, 
corresponding to a co-moving observer and k is the intrinsic spatial curvature. 

For the above FRW model of the universe if we make a change of the radial coordinate : r —> R where 
R = ar is known as area radius then the above standard FRW metric become Pinleve-Gustrand-like metric [3D] 
as follows 

ds 2 = --^JU 2 - ^^dtdR + T ^ + R 2 dn 2 (2) 

a 2 a 2 a 2 

Here Ra = — , 1 is the radius of the apparent horizon. 

v ff2+ ^ 

In the present case for the metric ^ the Kodama vector and the corresponding energy are 



k»= I 0, 0, 



and 



.,_ 1 1 kR 2 dS 



(3) 



Thus — -J ^ kR .j is the energy of the particle as measured by an observer with the Kodama vector. 
The differential equation for the radial null geodesic (i.e., ds 2 = = dQ 2 ,) has the form : 




™ + (4) 
R A 

where +/— sign are associated with outgoing /ingoing null geodesic, with the assumption that 't' increases 
towards future. As we are interested in the imaginary part of the action corresponding to the tunnelling process 
through a barrier (the classically forbidden region) so according to Parikh-Wilczek [T9] 

ImS = Im p R dR = Im / dp' R dR = Im / ^dR (5) 

JR in JR in JO JR irl JO R 
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where in the last step we have used the Hamiltonian equation 

^ _dH _ dH 
dpR dp R 

and here pr is the radial momentum i?i„, Rout are positions very close to horizon with Ri n the initial position 
and Rout, a classical turning point. Using the value of R from equation (j4|) into ([5]) we get 



r R, 

ImS = Im I dR 



dH' 

R in J R 



R ^ dR 



-U) 



= Im 

R ° ut dR 
Ri ' 



- v /r +^{v w + 1 -f - HR } 

= itRau (6) 
where integration over the Hamiltonian 'if' gives the energy of the particle measured by an 

observer with Kodama vector. 

Now comparing the tunnelling probability 

r^expj-^lmsj (7) 

with the Boltzmann factor exp{ — we have the temperature associated with the apparent horizon of the 
FRW universe as 

This is the semi-classical Hawking-like temperature of the FRW universe for tunnelling of massless particles 
across the apparent horizon. 

3 Hamilton-Jacobi method : Quantum prescription 

This section deals with tunnelling of massless particle beyond semiclassical approximation by Hamilton-Jacobi 
(HJ) method. 

We start with the Klein-Gordan(KG) equation for a scalar field <j> describing a massless scalar particle of 
the form 

h 2 



=a„ (s M V=^) = (9) 

The explicit form of the KG equation in the background of the above FRW metric (pj is given by 

dt 2 a? dr 2 dt a 2 dr 

It should be mentioned here that due to spherical symmetry of the FRW space time and consideration of the 
radial trajectories only we have considered (t — r)-sector in the space time given by equation ([1]), i.e., 2D 
hyperplane (t, r) . Now substituting the standard ansatz for the semiclassical wave function 

(»■■ /) <'xi>^s(r, t)\ run 



3 



into the wave equation (|10p. the differential equation for the action S becomes 



as 

~dt 



1 - kr 2 



dS_ 

dr 



ifi 



d 2 S 
dr 2 



H 



dS_ 
~dt 



kr dS 
a 2 dr 







(12) 



As a first step to solve this partial differential cquation(p.d.e.) we expand the action in powers of Planck 
constant % as 

S(r, t) = S {r, t) + X k h k S k (r, t) (13) 

with k, a positive integer. Here terms of the order of Planck's constant and its higher powers are considered 
as quantum corrections over the semiclassical action So. Now substituting this ansatz for S in the differential 
equation (|12l) and equating different powers of % on both sides we obtain the following set of partial differential 
equations : 




.and so on. 



= 



Apparently, different order p.d.c.s are very complicated but fortunately there will be lot of simplifications 
if in the p.d.e. corresponding to order h k , all previous p.d.e.s are used and finally we obtain identical p.d.e.s 
namely 



0S k \ 2 (l-kr 2 \(dS k ^ 2 



dt 



dr 







fc = 0,l,2,. 



(17) 



We see that different order quantum corrections satisfy identical differential equations as the semiclassical 
action So, so the correction terms are not independent, rather proportional to So (i.e., S k oc So for all k). 
To determine these proportionality constants we shall use dimension analysis. Since 5*0 has the dimension h so 
the proportionality constant d k for S k has the dimension h~ k . But in standard units namely G — c = ks = 1, 
the Planck's constant fi is of the order M 2 (M p =Planck's mass) and hence d k has the dimension M~ 2k , where 
M is identified as the mass of the universe. Thus the scries expansion (f]~3|) can be written in terms of So as 



S(r, t) = S (r, t) 



1 + SfeQffe 



h 

M 2 



(18) 



with a k as dimensionless constant parameters. 

Now, a complete solution for S requires solution of So satisfying the p.d.e (fT4)l . Since the metric (pj is 
non-static so there is no time like Killing vector in the dynamical FRW space-time. However, Kodama vector 
[55] has similar role in FRW spacetime as the timelike killing vector does in the stationary black hole space 
time. Inside the apparant horizon Kodama vector is a time like vector and there is a conserved energy of a 
particle moving in the time like killing vector in the stationary black hole space time. 

The Kodama vector for the FRW metric is given by 



K» = (V 1 - kr 2 , - Hr\/l - kr 2 , 0, o) 
with associated conserved energy of the particle 



UJ 



- \J\ — kr 2 — — - + Hr\J 1 — kr 2 -—^ 
at or 



:dSo 



Solving for ^f- and ^r 2 - using equations and (|2T)]) we get 



5-0 = 



u>dt 



{vT— kr 2 — iJar} 



=F aui 



dr 



y/l — kr 2 { vT— kr 2 — Har} 



(19) 



(20) 



(21) 
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where — /+ sign corresponds to ingoing/outgoing scalar particle for which the wave functions have the expres- 
sions 



fan 

and 

4>ir, 



exp 



--< 1 + ^a k [ 



idt 



exp 



{Vl - kr 2 - Har} 



L)dt 

{Vl — kr 2 — Har} 



dr 



Vl — kr 2 { vl — kr 2 — Har} 



dr 



Vl - kr 2 {Vl - kr 2 - Har} 



(22) 



(23) 



Now, across the horizon the metric coefficients in the (r, t) sector alter their sign, so the above time integration 
may have imaginary part and have contribution to the probabilities for the ingoing and the outgoing particles. 
As a result the probabilities are given by 



p. 

J in ' 

and 

Pout 



\<pi 



exp 



1 + Sfcafc 



h 

M 2 



Im 



u)dt 



{ Vl — kr 2 — Har} 



Im aiu 



dr 



exp 



- { 1 + S fc a fc ( 



Im 



jdt 



{Vl — kr 2 — Har} 



Im auj 



Vl - kr 2 {Vl - kr 2 - Har} 
(24) 

dr 



\fl — kr 2 {Vl — kr 2 — Har} 
(25) 

But in the classical limit (h — > 0) the outgoing probability has to be unity as then there will be no absorber 
and everything will go out [55] and hence from (|2"5]1 we get 



Im 



jdt 



{Vl - kr 2 - Har} 
So Pi n will have the simplified form 



Imauj 



dr 



Vl — kr 2 {Vl — kr 2 — Har} 



(26) 



p 



exp 



exp 



4aw 
2w J 



1 + S fc a fc 



1 + Sfeafe 



h 

M 2 

h 

M 2 



Im 



ttR a 



dr 



Vl — kr 2 {Vl — kr 2 — Har} 



(27) 



Hence from the principle of 'detailed balance' , i.e., 

P out = exp 

we have 



"n 



p, 



1 + ^a k [ ^ 



1 



2irR A 



= { 1 + S fe a fe ( 



T, 



(28) 



This is the horizon temperature of the FRW model of the universe after quantum corrections are taken into 
account. The quantum corrections to this modified temperature is very smilar to that of Banerjee et. al. |26j 
for general static black holes. Also the modified temperature has arbitrariness due to choice of the parameters 
afe in the expansion. Banerjee et. al. has shown for static black holes that different choices of o^'s will lead 
to different physical interpretation. For future work, we shall attempt to study particles with no-zero mass and 
examine whether quantum corrections as well as Hawking like temperature depend on the mass term. 
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4 Hamilton-Jacobi method in painleve type coordinate system for 
FRW universe 

As in the previous section The KG equation for a scalar field <f> describing a massless scalar particle has the 
form 





h 2 Q0 



i.e., 



where the box operator is evaluated in the back ground of the metric ([2J in the (t, r) part of the space-time as 
before. The explicit form of the KG equation is 



C 



where c 



d 2 ib , d 2 ib 
_j_ 

dt 2 dr 2 



2B 



d 2 ip 
dtdr 



d (B 



B 



HR 



dR\/j, 



A 



d_ (C_ 
dt \fi 



dip 

at 



d (A 



d (B 



OR \ fi ) dt \ ji 



1 = ° M 



and // = AC + B 2 . 



l— ! HP l- 

If we substitute the standard ansatz for semi classical wave function namely 



ij)(R, t) = exp t) 

then the action S will satisfy the following differential equation 

dS dS 



(30) 



C 



ifi 



A 



d 2 S 
dR 2 



C 



d 2 S 
dt 2 



2B 



d 2 S 
dRdt 



d 2 s 


-A 


d 2 S 


dt 2 


dr 2 








— 

OR 




)- 



25- 



d (B 



dt \ fi 



dt OR 
dS_ 
dR 



d (B 



dR\\i 



d_ (C_ 

dt \ Jl 



dS_ 
~dt 



= 
( 31 ) 

In order to incorporate quantum corrections over the semiclassical action, S will expanded in powers of Planck 
constant h as 

S(R,t) = S (R,t) + T,h k S k {R,t) (32) 



where Sq is the semiclassical action and K is a positive integer. Substituting 
powers of h on both sides we have the following differential equations : 



dSo 
dR 



+ ih 



A 



h 1 : 

d 2 S 
dR 2 



dSodS 1 _ dSodS 1 
dt dt dR dR 



2B 



dR dt 

fdSodSi dS dS! 
\~dR~dT + ~df ~dR 



in pip and equating different 

(33) 



- B- 



.d 2 Sn 



dt 2 



-2B 



d 2 S 
dRdt 



A 



\dR\n 
and so on. 



a_ /b 

dt 



dSo 
dR 



{-(- 
\dR\n 



d_ (C_ 
dt \Ji 



dSo 
dt 



= 
(34) 



To solve the semi classical action So from equation we first write in terms of ^f- as 



dS 
dR 



HR± 



H 2 R 2 + (] 



R 2 

7^ 



dSn 



R 2 



dt 



(35) 



where +/— sign corresponds to ingoing/outgoing scalar particle. For the Painleve-Gulstrad type coordinates 
given in equation ([2} the energy of the scalar particle can be obtained using Kodama vector [3S] as 
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LO = —\ 1 — 



kR 2 dS Q 



(36) 



where 1 — is the Kodama vector for the metric ([2]) . Using the value of ^f- from (|36|) we obtain 

using ([35]) and hence 



S (R, t) 



Lodt 



1 



R 2 



1 - 



fl2 



1 



kR 2 



(37) 



Then the wave functions denoting ingoing and outgoing solution of the K.G. equation (|29| using (|30|) are of the 
form 



and 



V>m = exp 



ip ou t = exp 




(38) 



(39) 



As in the previous section due to tunnelling across the horizon the coordinate nature changes, i.e., the sign of 
the metric coefficients in the (R, £)-hyperplane are altered so the time integral has an imaginary part and hence 
it has a contribution to the probabilities [26, 37 . Thus both incoming and out going probabilities are given by 



Pin = \ipin\ = exp 




Lolm 



HR+Jl 



R? 



AR 



1 - 



kR 2 



(40) 



and 



Pout = \tpout\ = exp <( - 



Im 



\ 



Lodt 



— ujlm 



kl! 2 



HR+Jl -j£ 



R 2 



-.dR 



kR 2 



(41) 



Now some simplifications are possible using the fact that all outgoing particles certaining cross the horizon, i.e., 
Pout = 1. So from equation (|4*T1) 



Im 



Lodt 



= oolm 



1 - 



kR 2 



HR- Jl - 



R 2 



-AR = 



1 - 



kR 2 



(42) 



and so Pi n simplifies to 



Pin = \(j>in\ = exp { —Im 



HR+Jl 



R? 



-AR 



i. ii- 



= exp < -j7-KR A 

Then from the principle of " detailed balance" [22j [231 El we write 



(43) 



exp 



P out = exp 



LO 



(44) 
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Thus comparing (|43|) and (l44l) , the temperature of the horizon is given by 



Hence we have the same expansion for the temperature as we have obtained in previous sections. Thus choice 
of coordinate is not matter for the HJ method. Here higher order correction terms can be obtained in a 
straightforward manner and lead to identical results as in the previous section. 



5 Determination of Entropy : Area law 

I this section we shall examine the semiclassical Bekenstein-Hawking area law [351 E3] for the non-static FRW 
space-time at the apparent horizon. Also we shall calculate the corrections to the semiclassical entropy due to 
quantum effects of the Hawking like temperature. So we start with the thermodynamical law which express the 
energy conservation as 

dM = T h dS A (46) 

Here Th is the temperature of the horizon (with quantum corrections), Sa the entropy of the horizon and M is 
chosen as the Misner-Sharp gravitational mass, defined as [7] 



(47) 



on the horizon. 

Now if we identify M in the quantum corrections (|28[) for the temperature as the above Mishner-sharp mass 
then integrating (|46| we have 



S A 





f 4ttM 


I™- 

I -Lh j 


/ h 



1 + 



h 

M 2 



dM 



— 1- 47rai In M - 4tt^- + O (h 

M 2 

= + 47rai In A/ - 16tt^ + O (h 

R A 



h 

1 A 2 

2 An 



(48) 



In the above the first term on the r.h.s. is the semiclassical relation between area and entropy which is the usual 
Bekenstein-Hawking area law with a discrepancy of factor ~. The second term is the leading order quantum 
correction term and is the standard logarithmic correction term in black hole thermodynamics [HJ 2SJ [33 HZ] ■ 
The higher order correction terms are in the inverse powers of the area of the horizon. 



6 Summary 

In this work, we have studied Hawking like radiation from the FRW model of the universe using both the ap- 
proaches — the radial null geodesic method (tunnelling approach) and the HJ formalism. In both the methods 
we have obtained identical Hawking-like temperature at the semiclassical level. Also we have shown that choice 
of coordinate is not at all matter in HJ approach. The factor of 'two' problem in the tunnelling approach has 
been overcome considering Kodama vector instead of time-like vector and associated energy as the energy of the 
tunnelling particle. We have obtained quantum corrections to the Hawking-like temperature using HJ method 
and these corrections are similar to these for general static black holes. Also, quantum corrected entropy for- 
mula has been evaluated from the law of thermodynamics and it is found that usual entropy law with quantum 
corrections are same as in black hole thermodynamics with some discrepancies in the multiplicative factor. For 
future work we shall attempt to find interpretation of the parameters involved in the quantum corrections. It is 
worthy to generalize the tunnelling approach for non-zero mass particles(i.e., time-like geodesies) as well as to 
incorporate the quantum correction terms (back reaction effects). Finally, we shall try to resolve the ambiguity 
in the multiplicative factors for the entropy-area formula. 
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